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We investigate how well the relic density of dark matter can be predicted in 
mSUGRA. We determine the parameters to which the relic density is most sensi- 
tive and quantify the collider accuracy needed to match the accuracy of WMAP 
and PLANCK. 



One of the very attractive features of the minimal supersymmetric stan- 
dard model is that it provides a natural candidate for cold dark matter. 
Cosmology, and in particular the measurement of the relic density of dark 
matter then put some strong constraints on the supersymmetric model. For 
example, WMAP [1] which at 2a constrains the relic density in the range 
.094 < Vlh? < .129 points to very specific scenarios [2]. Within the context 
of mSUGRA models, the favoured scenarios are slepton coannihilation, heavy 
Higgs annihilation and Higgsino LSP. The latter occurs in the focus point re- 
gion with heavy scalars. Predictions for the relic density depend sensitively 
on the details of the supersymmetric models but also on some assumptions on 
the cosmological model. Assuming the LHC will discover supersymmetry, the 
question then is whether one can make precise predictions on the value of the 
relic density of dark matter, hence confront the cosmological model. 

In the following we investigate to which extent precision measurements at 
a linear collider can match the accuracy on the relic density of dark matter as 
measured my WMAP (10% level) or as expected from PLANCK (2% level). 
To estimate this we will concentrate on a few observables, the ones that are rel- 
evant in the rather fine-tuned mSUGRA scenarios allowed by WMAP without 
assuming mSUGRA when analysing data [3]. The relic density is computed 
with micrOMEGAsl . 3 [4] and the evaluation of the supersymmetric spectrum 
relies on Sof tSUSYl.8.7 [5]. 
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1 Coannihilation 

In mSUGRA at small M there exists a region with almost degenerate fi — Xi- 
The contribution of coannihilation channels to the effective annihilation cross 
section, in particular Xi^i — * T 7 or fi^i — * tt j is essential in bringing the 
relic density in the desired range. In computing the effective annihilation 
cross section, coannihilation processes are suppressed by a Boltzmann factor 
oc exp~ AM / Tf where AM is the mass difference between the NLSP and the 
LSP and Tf the decoupling temperature. One expects Wi 2 to be very sensitive 
to this mass difference. Following Ref. [6], we take a slope "SI" in parameter 
space with \i > 0, tan/3 = 10, A — and 

M = 20.3333 + 0.134M 1/2 + 5.0667 x l(T 5 (M 1/2 ) 2 (1) 

where masses are in GeV. Along this slope (Mw 2 = 350 — 950 GeV), the relic 
density is in rough agreement with the WMAP range. Fig. 1 shows that the 
mass difference between fi and Xi varies from 12 GeV to 300 MeV as one 
increases the f\ mass. The mass difference between the LSP and the sleptons 
of the first two generations can be much larger. To estimate the sensitivity of 
flh 2 on the mass difference we have, for each point on the slope, varied the mass 
of the fi while keeping all other parameters fixed. We find that for the staus 
accessible at LC500, the mass difference must be measured to 0.9(0.2) GeV 
to match the WMAP (PLANCK) accuracy as displayed in Fig. la. Recently 
two groups have studied the detectability of nearly degenerate fi — Xi as wei l 
as the feasability of a precise measurement of the fi mass at linear colliders 
[7] . The results indicate that the level of accuracy needed to match PLANCK 
might just be reached. 

2 Large tan/3 

In the Higgs funnel region, rapid annihilation of neutralinos proceed via the 
s-channel exchange of a heavy pseudoscalar Higgs. This region is found at 
large values of tan/3 and is accessible to LC500 only for low values of M 1 / 2 - 
Within mSUGRA the relic density is sensitive to both mb(m,b) as well as the 
top quark mass since these strongly affect the value of the pseudoscalar mass. 
One can reduce significantly this theoretical uncertainty by measuring directly 
some of the physical parameters. The crucial parameters in this scenario are 
the masses ttl^o^Ma- As an example we consider a mSUGRA model with 
tan/3 = 50, Mi /2 = 380GeV, M = 399GeV, A = > 0, and then vary 
either m^o (through Mi) or Ma while keeping all other parameters fixed. We 
find that to match WMAP (PLANCK) accuracy on flh 2 , M A and m^o need 
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to be measured with a precision of 1%(0.2%). Note that the LHC can measure 
precisely Ma provided the A — ► channel can be used [8]. 

3 Higgsino LSP 

In mSUGRA, a LSP with a significant Higgsino component can be found at 
high values of Mo near the boundary of viable electroweak symmetry breaking 
where the value of fi drops rapidly. This is the focus point region. The en- 
hanced coupling of the Higgsino pairs to the Z and to the Higgses then makes 
for efficient annihilation into gauge bosons or ti pairs. Coannihilation channels 
with X2/X1 & l so contribute to the relic density. 

The position of the focus point region is extremely sensitive to the value 
of the running top quark mass that enter the spectrum calculator codes [3,9]. 
For m t — 175 GeV, the region where fj, « Ml and where the LSP has some 
Higgsino component is found generally for high values of tan (3. We will there- 
fore consider a model with tan/3 = 50, fx > 0, A = 0, and define the slope S2 
where the value of Slh 2 is in rough agreement with WMAP: 

m = 3019.85 + 2.6928M 1/2 - 1.01648 x 10~ 4 (M 1/2 ) 2 . (2) 

We have estimated the sensitivity of the relic density prediction to m t along S2 
in mSUGRA. Reaching WMAP precision requires knowing m t to 20 MeV. This 
is too demanding considering that at a linear colider one expects a precision 
of lOOMeV on m t when including theoretical uncertainties. 

The situation improves significantly if one does not rely on mSUGRA and 
uses the information that could be provided by colliders, the mt dependence 
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drops to the few percent level. Here the Higgsino fraction determined by n, Mi 
is expected to be the crucial parameter. Along the slope S2, we vary either /u, 
or Mi and recalculate the whole spectrum. We then estimate the accuracy a 
required on each parameter to get a 10% shift in flh 2 . As displayed in Fig. lb, 
a precision of 1 — 1.5% is required, while matching PLANCK accuracy requires 
roughly an accuracy of 0.2 — 0.3%. Since m^o and m^a are directly related to 
H and Mi a precise measurement of these two masses would allow to recon- 
struct the necessary parameters. The situation for tan j3 is more problematic 
as even the precision required for WMAP (10%) might be difficult to reach. 
Nevertheless a realistic study of the achievable accuracy in this scenario at a 
linear collider still needs to be performed. 

4 Conclusion 

The precise knowledge of the cross-sections necessary to make an accurate 
prediction of the relic density of dark matter rests on a precise knowledge of the 
physical parameters of the MSSM. Although the relic density calculation often 
involves a large number of processes, within the scenarios that are favoured by 
WMAP, only a few parameters need to be measured with very high precision. 
Furthermore, using the collider data rather than relying on some theoretical 
prejudice, considerably improves the precision of the prediction of Qh . 
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